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Abstract 

A new strategy to control the flow of surface plasmon polaritons at metallic surfaces is 
presented. It is based on the application of the concept of Transformation Optics to devise 
the optical parameters of the dielectric medium placed on top of the metal surface. We de- 
scribe the general methodology for the design of Transformation-Optical devices for surface 
plasmons and analyze, for proof-of-principle purposes, three representative examples with dif- 
ferent functionalities: a beam shifter, a cylindrical cloak and a ground-plane cloak. 
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Transformation Optics (TO) has been proposed^^ as a general technique to design complex 
electromagnetic (EM) media with unusual properties. In order to construct novel optical devices, 
one has to imagine a space with some distortion in a certain region and make a coordinate transfor- 
mation from a flat empty space. Due to the form-invariance of Maxwell equations under general 
coordinate transformations,!^ the properties of the distorted geometry can then be interpreted as 
that of a medium in the original flat empty space.l^ Thus, TO provides us with expressions for the 
dielectric permittivity tensor, e, and the magnetic permeability tensor, fl, that need to be imple- 
mented in order to obtain a medium with a designed functionality. The first practical realization 
of this idea was the construction of a two-dimensional (2D) cylindrical invisibility cloak for EM 
waves in the microwave regime. ^ This experiment took advantage of the recently developed field of 
Metamaterials. 1^ These are artificially structured materials made up of subwavelength constituents 
and designed to implement a prescribed response to EM fields. The cylindrical cloak experiment 
was followed by the construction of broad-band cloaks in the microwave*^ and optical*^^ regimes, 
all of them based on a different strategy, the so-called ground-plane cloak.^Hl Besides, a wide va- 
riety of applications other than cloaking has been recently presented, including beam shifters and 
beam dividers, field concentrators collimators ,^21 waveguide bends and comers,^ lenses with 
subwavelength resolutionl^^ or deep-subwavelength waveguides.'^ 

In Plasmonics, one of the main goals is to control the flow of light at a metal surface by means of 
the surface plasmon polaritons (SPPs) that decorate a metal-dielectric interface.^^^^ In the quest 
for creating photonic circuits based on SPPs, different waveguiding schemes have been tested 
during the last years (see, for example, the recent review articles™!). At a different level, 2D 
optical elements for SPPs like mirrors and beam splitters based on Bragg-scattering of SPPs with 
periodic arrays of scatterers have been also devised.!^ 

In this paper we present a different strategy to tackle this problem by showing how the TO 
framework can be applied to mold the propagation of SPPs at a metal surface. As TO is valid for 
any source of illumination due to the universal character of Maxwell equations, it is expected that 
TO recipes will be also operative for SPPs. These EM modes are surface waves travelling along 
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the interface between two media, a dielectric and a metal. Therefore, in order to design an optical 
device for SPPs, the expressions for the EM material parameters provided by the TO formalism 
should, in principle, be implemented both in the dielectric and the metal sides. This fact is a 
severe technical challenge, as TO generally requires highly anisotropic and inhomogeneous £ and 
jU. Moreover, the decay length of the SPP within the metal, i.e. the skin depth, is subwavelength, 
which means that a manipulation of e and /i at a nanometer scale within the metal would be needed. 
However, we will show in the following that a simplified version of the TO recipe in which only the 
dielectric side is manipulated leads to quasi-perfect functionalities. Two factors determine whether 
this simplification is operative or not: the wavelength of the SPP and the geometry of the device 
under consideration. In order to illustrate this point we discuss here three very different devices: a 
beam shifter, a cylindrical cloak and a ground-plane cloak for SPPs. 

As a first illustration of our methodology, we consider the 3D implementation for SPPs of 
a parallel beam shifter recently proposed by Rahm and co-workers for a 2D geometry.^ Let us 
discuss first this 2D case that was the first report of a finite embedded transformation, which is 
capable of transferring the modification of the EM fields to the wave that exits the TO medium. 
In the design of the 2D beam shifter, a rectangular region of sides d and / is considered and a 
transformation is performed from a 2D Cartesian grid into a grid that is tilted at an angle ^ [see 
inset of Figure 1(a)]. The amount of shift that the incident beam will experience is given by 
a = tan(^, which we will set to 1 = 45^) throughout this work for simplicity. The symmetry of 
the 2D problem for a EM plane wave implies that the two polarizations can be treated separately. 
For a wave polarized with the electric field perpendicular to the plane of the transformation, a TM 
plane wave, the relevant EM parameters obtained for the transformation medium according to TO 
are e^z = 1, fJ-xx = 1, Mxv = Myx = a and /iy^ = 1 + a^. 

To design a 3D beam shifter for a SPP propagating along the interface between a metal and a 
dielectric, we perform the embedded transformation in parallelepipeds of rectangular base (d x I). 
Two different transformation media are needed, one in the dielectric, with height h^, and one in 
the metal (height hm)- Notice that in order to operate over the entire SPP field, the height of the 
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rectangular slab in each medium must be greater than the corresponding decay length of the SPP 
in that medium (metal or dielectric). The expressions for £ and jX inside the dielectric slab in a 
cartesian basis are independent of the z coordinate and are: 



£to — f^TO 
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(1) 



where we have assumed that the dielectric half-space is vacuum without any loss of generality. On 
the other hand, the two tensors in the metallic slab should be: 



£ = £ro ■ A = Aro (2) 

where £„, is the electric permittivity of the metal. 

To test the TO recipe for a SPP beam shifter, we have carried out 3D full-wave simulations 
using the COMSOL Multiphysics software that is based on the finite element method (FEM). 



[figure] [1][]1 ^a) shows a 3D view of the behaviour of a SPP at A = 1.5 /im going through a beam 
shifter. The incident SPP, which displays a gaussian profile in the j-direction, propagates along the 
x-direction through a gold- vacuum interface and enters into the transformation medium, where it 
is shifted in the y-direction. Details on the geometrical parameters can be found in the caption of 



[figure] [1][]1 In this simulation, rectangular slabs with the prescribed £ and /i were introduced 



in both the dielectric and metal sides. A cross-section showing the behaviour of Hy at the metal 



surface (z = 0) is rendered in [figure][l][]l ^b). [figure][l][]l c) and (b) allow us to compare the 



quality of the shift whether we implement £ and A in the dielectric and metal sides (b) or only 
in the dielectric side (c). It is worth noticing that there is almost no difference between the field 



distributions in both cases. On the other hand, [figure] [1][]1 ^d) and (e) show the corresponding 



cross-sections at z = of the beam shifter for a SPP but now operating at A = 600 nm for the 
case of placing the transformation media in the dielectric and metal sides (d) and when only the 
dielectric side is manipulated (e). Although the effectiveness of the shifter at this wavelength is 
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affected by the absence of a transformation medium in the metal side, its functionality (a shift 
in the SPP beam) is still very good. The reason for this weak wavelength dependence is clear. 
The vacuum decay length of a SPP on a gold surface at 1.5 jim is 2.4 jim, much larger than the 
skin depth (23 nm) so the EM fields within the metal carry 0.01% of the total SPP energy, having 
practically no relevance compared to the EM fields in the dielectric. However, at A = 600 nm, the 
ratio between vacuum decay length of the SPP and skin depth is reduced by a factor of 10 and the 
EM fields inside the metal are then relatively more important for the propagation of the SPP. For 
this wavelength, 1% of the SPP energy resides in the metal. 

As a second example, we present a 3D cylindrical cloak for SPPs based on TO. Here we aim 
for suppressing the scattering of a SPP from an object placed on a surface and not to use plasmonic 
structures to achieve cloaking.l^^ Note that other strategies for SPP cloaking that do not rely on 
the TO framework have been recently proposed.'^ As stated above, the 2D invisibility cloak was 
the first application of TO, both in theoretical proposalsl^^^ and experiments .'2' The design of a 3D 
cylindrical cloak for SPPs is analogous to that of the beam shifter. The field distribution of the SPP 
extends itself in the dielectric and in the metal, so, in principle, we should again cloak the SPP-field 
in both regions of space. Because of the presence of the metal-dielectric interface, we apply the 2D 
cylindrical transformation as described by Pendry and co-workers ^ in planes that are parallel to the 
interface, building-up two concentric cylinders of radii a and b [see Figure 2(a)]. In each plane, we 
perform a transformation from a flat 2D euclidean space, where light propagates in straight hues, to 
a space with a hole of radius a surrounded by a compressed region of radius b. The transformation 
medium is the shell comprised between a and b and is expected to guide the incident SPP around 
the hole. The SPP-wave will emerge at the other side as if it had travelled through an empty space. 
The e and fl tensors obtained for the transformation medium in the dielectric half-space read as 
follows in a cylindrical basis: 
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^ '-^ ^ ^ 

ho = V^TO= (3) 

Cloaking the SPP field inside the metal would require a cylindrical transformation medium whose 
magnetic permeability is the same as in the dielectric side, /iro^ while its electric permittivity is 
simply £TO-£m- 

3D simulations of the cylindrical cloak for SPPs operating at A = 1.5 jim [panels (b) and (c)] 



and X = 600 nm [panels (d) and (e)] are shown in [figure] [2] []2 In these cases, the prescribed e and 
fl were implemented only in the dielectric, in a cylindrical shell surrounding a perfectly conducting 
cylinder, leaving the metal as a continuous sheet. We present for each wavelength the y-component 
of the scattered magnetic field evaluated at two planes parallel to the metal-dielectric interface, one 
at z = [panels (b) and (d)] and another at a height at which the SPP has decayed to a quarter of its 
intensity at the interface [panels (c) and (e)]. At A = 1.5 jim, the scattered magnetic field is nearly 
zero out of the cloaking region and we can conclude that the performance of the cloaking device 
for SPPs is excellent. However, for a SPP at a wavelength of 600 nm, only partial cloaking for a 
SPP is achieved. The reason for this is the geometry of the cylindrical cloak which, as opposed to 
the beam shifter, has a discontinuity in the metal-dielectric interface. In this case, the behaviour of 
the cloak is very sensitive to the field propagating inside the metal because it strongly scatters at 
this discontinuity. Simulations in which the metal is also manipulated according to the TO recipe 
were carried out (not shown here). In this case, the 3D cylindrical cloak for SPPs displays a perfect 
performance at the two wavelengths considered. Our results explain why a previous experiment 
on SPP cloaking at A = 500 nm in which the metamaterial was only placed in the dielectric side 
was not completely successful.!^ 

Finally, we consider a different strategy for cloaking,'I2based on the use of TO to design a cloak 
that mimics a ground-plane (flat metal surface). We will show that in this case the simplification 
of using the transformation media only in the dielectric side works perfectly. The main difference 
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between a ground-plane cloak and the rest of the TO-based devices that we have considered up 
to here is the geometry. A ground-plane cloak inherently needs to be placed on a metal surface, 
because it mimics a highly reflective surface instead of an empty space. Taking advantage of the 
presence of this metal-dielectric interface, we apply directly the TO recipe previously derived for 
a 2D ground-plane cloak illuminated by a plane wave'^to the case of SPP-illumination. 

The design of a 2D ground-plane cloakEl starts by considering a rectangular region of a flat 
space. The size of the rectangle isl xh and it sits on a metal surface. Suppose that we want to cloak 
a bump of maximum height ho and with a shape that follows the equation z{x) = ho cos^ ( jx) . Then 
we make a transformation from the rectangle to a region of the same shape except for its bottom 
boundary which is lifted following the surface of the bump. The simplest transformation consists 
of a transfinite mapping, which is not conformal.^^ This means that it does not preserve angles 
and hence the anisotropy is not minimized.'^ However, we will use it as a proof-of-concept of 
the performance of a ground-plane cloak for SPPs. If we place the cloak in the x — z plane (see 



[figure] [3][]3), the EM fields along the y-axis being translationally invariant, the transformation in 



a cartesian basis reads: 



x=^x',y = y',z = z +ho ^1 - cos^ (^jx' 



(4) 



The optical parameters, £ and jl, for the transformation medium are obtained using the standard 
TO techniques: 
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where A= fx^= ^ and f^^ = f^z^ + A^. The 2D ground-plane cloak for SPPs consists of a 
dielectric with the prescribed £ and fx placed on a metal surface with a cos^-shaped bump. A SPP 
propagating along the metal-dielectric interface and entering into the ground-plane cloak should 
be smoothly guided around the bump, [[figure] [3] []3[ a) shows the magnetic field distribution of 



a SPP at A = 600 nm travelling along a gold-vacuum interface and going through the cloak. Al- 
though only the dielectric medium is transformed in this 2D simulation, the SPP nicely follows the 
curvature of the bump at the two sides of the interface. 

A ground-plane cloak for SPPs can be also devised for hiding a 3D bump. In this case, the 
bump sits on thex — y plane and its surface is given by z{x,y) = /zqcos^ [jx) cos^ (fj)- As for the 
2D cloak, we make a transfinite transformation from a rectangular prism to a prism whose bottom 
boundary is lifted by the bump. In this case, the mathematical formulae for the transformation are: 



x = x' ,y =y' ,z = z +ho{ 1 - |- j cos^ ( ( "T"^' 



According to the TO rules, we obtain the e and fi for the transformation medium: 

( i A, ^ 
e = A = ^ 1 

^ fxz fyz fzz J 



(6) 



(7) 



where A = ^, f,^ = f^^ = ^ and = fj + fyJ + A\ We present in [figure] [3] []3 :b) a 



simulation that confirms the effectiveness of the 3D cloak for a SPP operating at A = 600 nm. The 
color map and the streamlines represent the density and the direction of the power flow of the SPP, 
respectively. The SPP is guided around the bump and continues travelling along the gold- vacuum 
interface. For comparison, [figure] [3] []3| ^c) shows a SPP scattering from the gold bump when the 
cloak is not present. We see that the fraction of the SPP that impinges directly into the bump is 
scattered upwards. This radiation into free space is totally suppressed when the cloak is present. 



as shown in |[figure][3][]3| ^b). 

Now we briefly discuss the practical realization of the TO-based devices for SPPs presented 
in this work. Our results have shown that the effectiveness of these devices is basically governed 
by the presence of a transformation medium in the dielectric side whereas modifying the metal 
side is not necessary in most of the cases. Then it is enough to realize anisotropic £ and fl within 
the dielectric, making use of a specific metamaterial. Different tools for overcoming the neces- 
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sity of anisotropic parameters have already been developed within the field of Metamaterials. For 
instance, in the reduced parameters approach, the tensors are simplified so that only some compo- 
nents need to be implemented. This allows the design of non-magnetic metamaterials operating in 
the optical regime, a fact that has been applied to the construction of the 2D cylindrical cloak. On 
the other hand, isotropic £ and fl are obtained when a conformal transformation^ is performed, so 
it is enough to implement an isotropic refractive index. Numerical techniques have been recently 
developed to compute quasi-conformal coordinate transformations, which minimize anisotropy to 
a point where it is negligible. In particular, the quasi-conformal mapping technique has proved 
useful for the practical realization of the 2D ground-plane cloak for plane wave illumination^ or for 
building-up a curved reflector that mimics a flat mirror for SPPs.l^ Similar ideas could be applied 
to the design of different devices, including the beam shifter analyzed in this work. Therefore, we 
are confident that the recipes that we have proposed could be practically implemented by properly 
structuring a dielectric layer placed on top of a metallic film at a length scale much smaller than 
the operating wavelength. 

In conclusion, we have demonstrated how the concept of Transformation Optics when applied 
to SPP propagation could offer the necessary tools to design new strategies for molding the flow 
of SPPs in structured metal surfaces. A simplification of the general TO recipes in which only the 
dielectric side is manipulated yields to almost perfect functionalities for SPP-devices. We hope 
that our results will encourage other experimental and theoretical groups to explore the exciting 
opportunities that the idea of Transformation Optics brings into the field of Plasmonics. 

Note: As we completed this manuscript, we accidentally found that another group^^also com- 
pleted a paper independently at the same time on a similar topic. 
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Figure 1: Normalized magnetic field pattern (color map) of a SPP travelling along a gold- vacuum 
interface and entering into a 3D beam shifter. Inset of the panel (a) shows the basic geometry 
of the 2D case. Panel (a) depicts a 3D view with power flow streamlines of a SPP operating at 
A = 1.5 jim. The prescribed £ and fl were implemented in rectangular parallelepipeds with in- 
plane dimensions 2.5/im x 10/im and thicknesses h^i = 6 jim and hm = 50 nm. Panel (b) renders 
a cross-section picture of panel (a) evaluated at the metal surface, z = 0. For comparison, panel 
(c) shows a SPP beam shifter where the slab with the prescribed e and fl has been introduced 
only in the dielectric side. Panels (d) and (e) render the operation of a beam shifter for a SPP of 
A = 600 nm. In this case the in-plane dimensions of the rectangular slabs are 1 /im x4/im and the 
thicknesses are = 600 nm and hm = 50 nm. In panel (e), only the dielectric side is manipulated. 
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Figure 2: Scattered magnetic field for a 3D cylindrical cloak for SPPs. (a) A schematic picture 
showing the geometry of the cloak. Cloaking parameters are implemented only in the dielectric 
side, in the region between the inner radius a and the outer radius b = 2a. The operating wavelength 
is 1.5/im in panels (b) and (c) and the cylindrical cloak in this case is 6jUm high with a = 1.5/im. 
(b) Scattered Hy field evaluated at the metal surface, z = 0. (c) The corresponding field as obtained 
at z = 625 nm, a quarter of the vacuum decay length of the SPP into the dielectric. In panels (d) 
and (e) the wavelength is 600 nm and the cloak is 600 nm high with a = 500 nm. (d) Scattered Hy 
field obtained at z = plane, (d) The same quantity but now evaluated at z = 70 nm, a quarter of 
the vacuum decay length of the SPP for A = 600 nm. 
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Figure 3: (a) Normalized magnetic field pattern of a 2D ground-plane cloak for a SPP travelling 
in the x-direction at 600 nm. The cloaked object is a cos^-shaped gold bump, 1250 nm long and 
with a height of ho = 200 nm. The height of the cloak is h = 750 nm. (b) Power flow distribution 
(color map) and streamlines of a 3D ground-plane cloak placed over a 3D cos^-shaped gold bump 
of dimensions 1250 nm x 1250 nm x 200 nm. (c) Power flow for the 3D gold scatterer without 
the cloak. 
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